Universal Scaling in Mixing Correlated Growth with Randomness 



in 

O 

o 

(N 

Cu- 
te 

vo 
(N 



c3 



-a 

o 
o 



> 

oo 

ON 

o 
in 
o 

-I— > 
c3 



i 

C 

o 
o 



X 



A. Kolakowska, M. A. Novotny, and P. S. Verma 
Department of Physics and Astronomy, and the ERC Center for Computational Sciences, 
P.O. Box 5167, Mississippi State, MS 39762-5167 
(Dated: February 2, 2008) 

We study two-component growth that mixes random deposition (RD) with a correlated growth 
process that occurs with probability p. We find that these composite systems are in the universality 
class of the correlated growth process. For RD blends with either Edwards- Wilkinson or Kardar- 
Parisi-Zhang processes, we identify a nonuniversal parameter in the universal scaling in p. 
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Many properties of complex systems can be un- 
covered by statistical analysis of some representative 
nonequilibrium interfaces. Mainstream studies of surface 
growth and interface roughening focus on one-component 
growth, or homoepitaxy, and large-scale properties. On 
the microscopic level, nonequilibrium interfaces have 
been studied in a variety of discrete simulation models 
such as ballistic deposition (BD), Eden or solid-on-solid 
models. While one-component growths are well under- 
stood the same can not be claimed about composite 
systems, even as simple as binary growth in one spatial 
dimension. Several mixed-growth models studied in the 
recent decade 0, reveal new and nontrivial properties. 
But the theory behind these is in the initial stages. 

In this first systematic study of two-component 
growth, we examine a system whose dynamics is gov- 
erned by two simultaneously present processes: one is 
a process that builds up correlations (a pure-correlated 
growth) and the other process is totally uncorrelated, i.e., 
random deposition (RD). The pure-correlated growth oc- 
curs with probability p. Questions that we address here 
concern the universality of such composite systems. As 
we shall show, the presence of randomness slows down the 
dynamics of the correlation processes. Nevertheless, the 
universality class of the combined processes is the same as 
the universality class of a correlation process. This is an 
outcome of scaling in p. One consequence of this observa- 
tion is a magnifying-glass effect that RD-blending has on 
the time-evolution of the surface roughness. This effect 
can be useful in revealing hidden features of a correlated 
growth when designing simulation models. Intuitively, 
since RD carries no correlations of its own, it may be 
expected that its admixtures should not lead to a new 
universality class. Yet, demonstration of this is not so 
trivial since, as we shall make evident by the results of 
several simulations, some of the parameters involved in 
the universal scaling may be nonuniversal. Results pre- 
sented here for (1 + 1) dimensions can be easily extended 
to multidimensions. 

Consider aggregation models where particles fall onto 
a one-dimensional substrate of L sites, where they may 
be accepted in accordance to a rule that generates cor- 
relations among the sites. This pure-correlated growth 



occurs with probability p and competes with RD growth 
that occurs with probability q = 1 — p. When a par- 
ticle is accepted at a site, the site increases its height 
by Ah. Roughness of the growing surface is mea- 
sured by the interface width w(t) at time t: (w 2 (t)) = 
(L^ 1 J2k=i [hk{t) — h(t)] 2 ), where hk(t) is the height at 
site k and h(t) is its mean over L sites (angular brackets 
denote the mean over N configurations). 

In a pure-correlated growth (p = 1), assuming elemen- 
tary linear and nonlinear models, the self-affined rough- 
ness obeys the Family- Vicsek (FV) scaling 



\t) = L 2a F(t/L z ), 



(1) 



where F(y) gives two evolution limits: F(y) ~ y 2a / z if 
y <C 1 (growth); and, F(y) ~ const if y ^> 1 (saturation). 
The cross-over time t x from growth to saturation is given 
by the dynamic exponent z, t x ~ L z (Fig. 1). At sat- 
uration the width does not depend on time, w 2 ~ L 2a , 
where a is the roughness exponent. During the growth 
w 2 (t) ~ i 2/3 , where [3 — a/z. Exponents z, a and (3 
are universal. This means, two different simulation mod- 
els of two different correlation mechanisms will generate 
the same type of scaling, with consistent values of ex- 
ponents, provided these mechanisms represent the same 
type of correlation process, i.e. , belong to one universality 
class. Dynamics of the buildup of correlations and dy- 
namical scaling are described within a continuum model 
by a stochastic growth equation. One example is the 
Kardar-Parisi- Zhang (KPZ) equation 



h t = v(t) + f h xx + {X /2)h 2 x + rj(x,t), 



(2) 



where h — h(x,t) is the height field (subscripts denote 
partial derivatives; x is the coordinate along the sub- 
strate) , v is the mean interface velocity, and r\ is the white 
noise (vq and Ao are coefficients). In the KPZ universal- 
ity class, governed by Eq. (J2J), a + z — 2 and a = 1/2. 
When Ao = 0, Eq. (J2J becomes the Edwards- Wilkinson 
(EW) equation |6(, defining the EW universality with 
2a + 1 = z and a = 1/2. When v a = A = 0, Eq. © 
describes uncorrelated processes of RD universality, char- 
acterized by/? = l/2,t x =oo, and the absence of scaling 
in L. For EW processes, Eq. (JJJ expresses the invariancc 
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FIG. 1: Model B for pure-correlated growth: (a) Time evo- 
lutions of the interface width (to marks the end of the initial 
nonscaling regime); (b) Scaling function for t > to- z, a and 
/3 are consistent with the EW universality. Here, N w 100. 
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FIG. 2: Scaled widths at saturation vs parameter 1/p 2 : (a) 
and (b) are for Models A and B, respectively (mix of RD 
with EW processes); (c) and (d) are for Models C and D, 
respectively (mix of RD with KPZ) . Reference lines have slope 
1. Data are scaled with the exponent values shown here. 



of the EW equation under the scaling |l( 

x-> Lx, h^> L a h , t-> L z t. (3) 

Similarly, for KPZ processes it expresses the invariance 
of the convective derivative in the Burger's equation. 

In simulations, t is the number of deposited monolay- 
ers. The first step is RD to a flat substrate. The system 
retains the memory of this initial condition for to steps, 
where to depends on the particulars of the model, i.e., to 
is a nonuniversal parameter. In this start-up regime w(t) 
does not scale 0; scaling occurs only for t > t (Fig. 1). 

In deriving the scaling hypotheses, we are guided by 
the following four models. Model A: for p = 1 is RD 
with surface relaxation where Ah = 1 Q, , known to be 
in the EW universality class (for p < 1, studied in 
In Model B, being introduced here, Ah is sampled from 
a uniform distribution of unit mean and the substrate is 
sampled sequentially at each t. When p = 1: particles 
that fall on the local interface minima are always ac- 
cepted; particles that fall on local maxima slide down to 
either of the neighboring sites with probability 1/2; and, 
particles that fall on local slopes slide down to nearest- 
neighbor sites. Model B for p = 1 is in the EW class 
(Fig. 1). It simulates, e.g., deposition of a sticky non- 
granular material of variable droplet size. Model C: for 
p = 1 is BD with Ah = 1, known to be of KPZ univer- 
sality 1] (for p < 1, studied in 0, In Model D, 
Ah is sampled from a Poisson distribution of unit mean, 
and each monolayer is obtained by sequential sampling. 
When p = 1 in Model D, particles are deposited only to 
local surface minima. This case is in the KPZ universality 



class [7|. Model D simulates, e.g., conservative updates 
in a system of asynchronous processors 0, El • We stress 
that, although in one universality class for p = 1, Models 
C and D are essentially different simulations (as is the 
pair A and B). 

In all models, evolutions w 2 (t) form two-parameter 
families of curves (L and p being parameters) that for 
any p £ (0; 1] look like those in Fig. 1 but with to = 
to(p) > io(l) 5 t x = t x (p) > t x (l), and at saturation 
w 2 {p) > w 2 (l). The curves saturate due to only one com- 
ponent, the pure-correlated deposition, since the other 
component, RD, introduces no correlations. At satura- 
tion, the observed lateral correlation length is £,\\(p) ~ L 

1/2 

and tx(p) ~ L z ; thus, £\\(p) ~ i x (p); and, the widths 
scale in L as w 2 (p) ~ L 2a . Plots of the scaled widths 
(w 2 (p))/L 2a (Fig. 2) show that they generally scale in 
p as w 2 ~ L 2a /p 2S , where S is some parameter. Is S a 
universal exponent? Models A and B (Figs. 2a-b) may 
suggest a universal value S — 1 for the RD-EW mix. But 
Models C and D (Figs. 2c-d) show that S D w 25 C - Ac- 
cordingly, S is not universal because for the RD-KPZ mix 
its value is clearly related to the technicalities of these 
models. In the RD-EW case there is no reason to believe 
that 5 = 1 is not accidental. Scale invariance of the EW 
equation is not sufficient to furnish S. 

Since t x (p) ~ L z and w 2 ~ L 2a /p 2S for any p ^ 0, the 
roughness must scale as w 2 (t)/w 2 ~ F(t(p)/ L z ). This 
scaling in L collapses all curves w 2 (t) to one-parameter 
families (p being the only parameter now) presented in 
Figs. 3a, 4a, 5a and 6a. As RD components do not build 
correlations, this collapse is obtained with the scaling 
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FIG. 3: Scaling for Model A: (a) in L; (b) in p of the data 
in Fig. (a), a and /3 are consistent with the EW universality 
class. Data labels are common for both figures. N ~ 100. 



FIG. 4: Scaling for Model B: (a) in L; (b) in p of the data in 
Fig. (a). As in Model A, a and (5 indicate the EW universality 
class. Data labels are common for both figures. TV ~ 100. 



laws from the corresponding universality classes of pro- 
cesses that build up correlations. Explicitly, z = 2a + 1 
and z = 2 — a for blending RD with EW and KPZ pro- 
cesses, respectively. To further collapse the data inp, i.e., 
to find t(p) in the argument of function F, we analyze the 
invariance of the corresponding continuum equations un- 
der simultaneous affine transformations: 

x-> Lx, h~> hL a /g, t->tL z /f , (4) 

assuming g and / being arbitrary suitable functions of p. 
Scaling is the superposition of scaling J3J with 

x' = x, h^ti = h/g(p) = t/f(p) . (5) 

Invariance analysis under the component scaling Q leads 
to Eq. JIJ and signature-scaling laws of KPZ and EW 
processes. This justifies the data collapse in L. The 
component scaling (J5J transforms Eq. to: h' t , = v' + 
v'(p)ti x , x , + (X'(p)/2)h'i,+r,'(x',t'), where h' t , = h t f/g, 
v' = vf/g, h' x , x , = h xx /g, and rj'(x',t') = yffr){x,t). Its 
invariance under © implies: 

M = 9 2 (P) (6) 

y'ip) = "of(p) (7) 
X'(p) = \ g(p)f(p). (8) 

From scaling at saturation we obtained g(p) = p s . Thus, 
the continuum equation for the RD-KPZ mix is 

ht = v(t) + wp 25 h xx + (X /2)p 3S h 2 x + n{x, t). (9) 

In the limits p — ► 1 and p — > Eq. © describes the 
dynamics of pure processes, i.e., the KPZ-type and RD, 



respectively. Similarly, the invariance of the EW equa- 
tion under scaling (JSJ gives Eqs. This leads to 
the continuum equation for the RD-EW mix: 

h t ^v(t) + ^p 25 h xx + r](x,t). (10) 

The inverse of the scaling Q is the desired contrac- 
tion that gives the full data collapse described by the FV 
function. The inverse of the scaling JSJ) alone (a; — ► x, 
h — > p s h, t — * p 2S t) gives t(p) in the argument of 
F(y) : V = p 2S t/L z . Finally, the FV scaling for any two- 
component process, one of which is RD, is 

= ^ iy) - (id 

where a and z are universal exponents of the component 
process that builds up correlations, and S is nonuniversal. 
This result is illustrated in Figs. 3b, 4b, 5b and 6b. 

Our results, Eqs. (|9"|)- (|llfl . show that mixing RD with a 
correlated growth preserves the universality of the corre- 
lated growth. Physical justification is in the uncorrelated 
nature of RD. As can be seen from Eqs. l(§)l- i|10|) . RD 
blending reduces the values of coefficients v and A rela- 
tive to the original noise strength. In other words, the net 
outcome is a noisier dynamics. The analysis presented 
here by the examples of EW and KPZ processes in (1 + 1) 
dimensions is easily extended to other growth processes 
in (1 + n) dimensions. It is enough to notice that Eq. © 
is generally valid when scaling (0) applies to growth equa- 
tions of the type ht(x,t) — (operator)Zi + rj{x,t), where 
x is n dimensional. Hence the conclusion: if a correlated 
growth belongs to a given universality class, its mix with 
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FIG. 5: Scaling for Model C. (a) in L; (b) in p of the data in 
Fig. (a), a and (3 are consistent with the KPZ universality. 
Data labels are common for both figures. Here, N « 100. 



FIG. 6: Scaling for Model D: (a) in L; (b) in p of the data 
in Fig. (a), a and j3 indicate the KPZ universality, but Sd ~ 
<5c/2. Data labels are common for both figures. N ~ 1000. 



RD will remain in the same class. The only effects of 
the RD admixture are the simultaneous dilatations of 
the fundamental time and height scales in accordance 
with scaling (J5J (and Eq. J§}). The net consequence of 
these is a slowdown in the dynamics of buildingup the 
correlations, reflected in the change of the lateral cor- 
relation length £||(p) - t(j>) x l z = £|| (l)/y/f(p)- In a 
sense, RD blending is like applying a magnifying glass 
to the evolution curves w(t): the smaller the p the bet- 
ter the magnification. In particular, in a two-component 
growth that mixes RD with either EW or KPZ processes, 
these dilatations explicitly are h — > h/p s and t — > t/p 2S , 
where S is nonuniversal and reflects the particulars of 
the deposition. The stretching in time causes the initial 
nonscaling regime fo(l) in curves w(t) to be amplified as 
to(p) = to{l)/p 2S ■ One consequence of this amplification 
is a clear observation of the RD growth (with (3 = 1/2) 
for initial times t < to when the growth starts from a flat 
substrate (e.g., observed in Note, if p — > this 

initial phase becomes infinitely long as this is the limit of 
RD growth. In simulations, when p is known, by a pru- 
dent design of a model, magnifying effects of RD blending 
may prove advantageous in revealing hidden features of a 
correlated growth. However, in the laboratory, the pres- 
ence of randomness in the growth process will obscure a 
clear-cut observation of the expected scaling. 
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